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We evaluate the pion and isobar propagators in cold nuclear matter self consistently applying a 
. covariant form of the isobar-hole model. Migdal's vertex correction effects are considered systemat- 

' ically in the absence of phenomenological soft form factors. Saturated nuclear matter is modeled by 

I scalar and vector mean fields for the nucleon. It is shown that the short-range dressing of the nNA 

, vertex has a significant effect on the pion and isobar properties. Using realistic parameters sets we 

predict a downward shift of about 50 MeV for the A resonance at nuclear saturation density. The 
pionic soft modes are much less pronounced than in previous studies. 
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I. INTRODUCTION 



The theoretical awroaches for the nuclear pion dynamics [H, i, H, i, i, 0, 0, i, B [13, Hll, [13, [H, [li, [H, H, 
[TtI . [TsI . [igl . [20I . [21I . [23 . [23} followed so far can be put into different categories. All works acknowledge and consider 
^ ' the important role of short-range correlation effects. However, there is no common consensus about their absolute 
^ I strength. The latter depends decisively on the subtle details of the considered approach. For instance the non- 
relativistic computations Til 13l. [l6t obtain contrasted results for the pion properties in cold nuclear matter. With 



few exceptions @, [T^, Il8l.fl9j non-relativistic many-body techniques are applied. Also works that incorporate the 
" ^ ' feedback effect of a dressed pion propag ator^ that depends sensitively on the isobar propagator itself, on the isobar 
' self energy are in the minority [sl. [T^. 1T4. [Tol [20l. [2ll. [23t . It has been found that self-consistency is a crucial effect for 
, the nuclear TriVA dynamics. Moreover, the in-medium isobar propagator should be used in the computation of the 
' isobar- hole contribution building up the short range correlation effects [l^ [l^ [l^l • In the early works that addressed 
] self consistency issues jl^ IT^ a quite soft phenomenological form factor was used. This implies a strong and artificial 
^\ • suppression of pionic soft modes with large momentum that dominate the isobar width [l9|. The use of such soft 
' form factors explains why in (l^ . [l^ quite conventional isobar properties [28j were obtained without the inclusion of 
00 , vertex correction effects in the isobar self energy. The use of soft form factors suppresses the in-medium mass and 
■ width shifts of the isobar significantly. Noteworthy is the work [Tsj , in which isobar properties were computed without 
^ , relying on soft form factors. The important role of a hard factors in the nuclear TriVA dynamics was pointed out in 
. One may conclude that a description of isobar properties [H, [13, [l^, [13, [Ml, [2^ that relies on soft form factors 
should not be considered microscopic unless one includes a strong density, energy and momentum dependence in the 
form factor. 

A further possibly important aspect is the splitting of the isobar modes in nuclear matter [3, fisl . [l^ . An isobar 
moving through nuclear matter manifests itself in terms of longitudinal and transverse modes described by distinct 
spectral functions. The splitting of the two modes was found to be small in [3, [ll]. In contrast, in [l3[ sizeable 
effects were found depending, however, on the precise structure of the form factors used. Notwithstanding, further 
clarification on the form of the isobar self ene rgy in nuclear matter is needed. This is of particular relevance for 
instance in applications to heavy-ion reactions [23| . 

Recently it was demonstrated [18| that a covariant form of the isobar-hole model differs significantly from non- 
relativistic versions thereof [3, [3,H]. Relativistic corrections are not important everywhere in phase space. As a 
striking example recall the behavior of the nucleon-hole contribution to the pion self energy. A proper relativistic 
treatment leads to a result proportional to uP' — q^ . with the pion energy and momentum uj and q respectively [gI. [sI. [18| . 
In contrast, a non-relativistic evaluation provides a factor (p only Obviously, the non-relativistic expression is 
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justified only in a small subspace of phase space. Paying contribute to this observation various prescriptions (see 
e.g. [l^l) were suggested in the literature. One may speculate that the incompatible treatment of such effects is an 
important source for conflicting sets of Migdal parameters used in the literature [1, [H, [13, [IS] • 

Though it should be possible to incorporate relativistic effects in a perturbative manner with possible partial sum- 
mations required, we argue that it is more economical to perform computations that are strictly covariant. Appl ying 
the projector techniques developed recently [isl. [2^. l25l. [26| it is straight forward to perform such calculations. In [l9| 
a first manifest covariant and self consistent computation of the pion self energy was presented. The incorporation of 
scalar and vector mean fields for the nuclcon was worked out recently in 1261 at hand of the nuclear antikaon dynamics. 

The purpose of this work is to extend the previous studies of two of us [l8l.[T9|. We compute the isobar self energy in 
a covariant and self consistent manner generalizing the covariant isobar- hole model of [3 ■ Vertex correction effects as 
well as the longitudinal and transverse isobar modes are treated consistently. Results will be presented for a range of 
parameters centered around a parameter set that was found to be compatible in [27t with the nuclear photoabsorption 
data [29j . The effect of various approximation is discussed and illustrated comprehensively. 



II. COVARIANT ISOBAR-HOLE MODEL 



We specify the isobar- hole model in its covariant form |17l . Il8| . The interaction of pions with nuclcons and isobars 
is modelled by the leading order vertices 



(1) 



where we use t} T, 



TiTj/i and /at = 0.988 and /a = 1.85 in this work. Short range correlation effects are 



modelled using the covariant forms of the Migdal interaction vertices as introduced in |17l . \li 



^Migdal = g'll^ [iJl5lp.fij) [lp-/5l''rij] 



(2) 



where it is understood that the local vertices are to be used at the Hartree level. The Fock contribution can be cast 
into the form of a Hartree contribution by a simple Fierz transformation. Therefore it only renormalizes the coupling 
strength in ^ and can be omitted here. The Lagrangian densities ([I] ^ are effective in the sense that we consider 
their coupling constants as functions of the nuclear density. This is justified since we do not incorporate the physics of 
higher lying baryon resonances nor further mesonic degrees of freedom like the vector mesons explicitly. Integrating 
out more massive degrees freedom leads to a density dependence of the coupling constants necessarily, which however, 
is expected to be quite smooth due to high-mass nature of the modes treated implicitly. 

Unfortunately, there is yet no set of Migdal parameters universally accepted. For instance, the computation [l^ 
used the universal values g'n = g'12 — g'22 — 0.60, based on a study of isobar properties. Universal values for the 
Migdal parameters were suggested first in Nakano et al. deduce the constraint ~ 0.585 together with 
g'12 = 0.191 -I- 0.051 g'22 insisting on the empirical quenching factor Q = 0.9 of the Gamow- Teller resonance [33- Their 
consideration assumes that the quenching results exclusively from a mixing of the nucleon-hole and the isobar-hole 
state. In our work the parameters g'^j are varied around the values g'n ~ 1.0, .922 = S12 — 0-4 obtained from a detailed 
analysis of the nuclear photo absorption data 

Our studies will be based on the in-medium nucleon propagator parameterized in terms of scalar and vector mean 
fields: 
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FIG. 1: Phenomenological isobar 
mass functions m'^'^{Vw^ ) that 
lead to the reproduction of the P33 
pion-nucleon scattering amphtude. 
The dashed (sohd) hne leads to an 
exact (approximate) reproduction 
of the amplitude. 



In the rest frame of the bulk with = (1,0) one recovers with (j4]) the standard result p = 2fc|,/(37r^). We assume 
isospin symmetric nuclear matter. 

The focus of our work is the study of the in-medium isobar propagator S^,j{w, u), the solution of Dyson's equation 

" ^ — niA + ie \ 3 3 3 ttja 

5^, {w, u) ^ 5^ {w ~^$u) + {w - S0 u) J:"f^{w, u) Sp,{w, u) , (5) 

where we allow for a vector mean field of the isobar. 

In nuclear matter the isobar self energy tensor, I]^i,(w, u), is a quite complicated object which involves the time-like 
4-vector characterizing the nuclear matter frame. In order to arrive at a reproduction of the P33 pion-nucleon 
partial-wave amplitude we allow for a phenomenological energy dependence in the free-space isobar mass. We write 



= (V^)-I]f, (6) 



where we introduce also a scalar mean field T,g for the isobar. At nuclear saturation density we found the values 
Tiy ~ —0.25 GeV and ~ —0.11 GeV to be consistent with the nuclear photo absorption data in an application 
of the present covariant and self consistent many-body approach (27| . Note that latter values are scheme dependent 
reflecting the particular in-medium processes taken into account explicitly. 

Making the assumption that the P33 amplitude is determined completely by the s-channcl exchange of the dressed 
isobar, for a given isobar self energy the mass function rn^^ can be expressed directly in terms of the empirical P33 
phase shift. Based on the self energy to be specified in section 5 we arrive at the mass function shown in Fig. 1 with a 
dashed line. The sizeable variation of m^'^ on VuP reflects contributions to the P33 amplitude that are characterized 
by left-hand branch points. The amplitude receives, besides the s-channel isobar exchange, additional contributions 
like from the nucleon u-channel process. Since the latter contribution will be considered being, implied by ([1]), it is 
not consistent to proceed with the dashed mass function of Fig. 1 . A fully consistent approach would require at least 
the unitarization of the sum of s-channcl isobar and u-channcl nucleon exchange processes. This is, however, beyond 
the scope of the present work. In order to correct for the presence of the u-channel nucleon exchange we determine 
the phenomenological mass function rri^'^{Vw^) in the following way: the s-channel isobar contribution is adjusted 
to reproduce the imaginary part of the P33 partial wave amplitude in the vicinity of the isobar peak. Away from 
the resonance the mass function is kept constant. The result is shown in Fig. 1 with a solid line. As compared to 
the dashed line, which reproduces the P33 amplitude exactly, the solid line shows a much reduced variation. This 
is welcome since the smoother the phenomenological mass function the smaller are the uncertainties implied by the 
ansatz ([5]). 

The quality of our prescription is illustrated in Fig. 2, where the empirical P33 partial wave amplitude in the 
convention of [l^ is confronted with the phenomenological amplitude. Real and imaginary parts agree well in the 
resonance region. Significant deviations are noted close to threshold only, where we expect a strong energy dependence 
from the u-channel contributions. This is confirmed by the additional solid line of Fig. 2 which shows the contribution 
of the u-channcl nucleon exchange process. Close to threshold it is largest almost making up the difference of the 
empirical and phenomenological amplitude. 
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FIG. 2: Comparison of the empir- 
ical P33 amplitude from [sj] with 
the phenomenological one as im- 
plied by ((61 with the solid mass 
function of Fig. 2. The addi- 
tional solid line is the contribution 
of the u-channel nucleon exchange 
process. 



III. PION SELF ENERGY 



In this section we evaluate the pion self energy as implied by the interaction ([T]) for a given isobar propagator 
S'^[/(u;,u). The latter will be specified in subsequent sections. The central objects to compute are the nucleon- and 



isobar- hole loop tensors, 11^^'* ■'(q, u) and Il\^'^'{q, u), which we define by 



{Ah) 



4 fl f d^P 



3 m? 



= 2 



■itr AS{p, u) S^u{p + q,u) + -q^) . 



+ ^'^S{p,u)-f5-ffj,AS{p + q,u)j5-/^ 1 + (9m ^ -9m) 



(7) 



with the isobar propagator, S^,^{w, u) of ([5]), and the in-mcdium part of the nucleon propagator, A.S{p, u) as specified 
in dSl). 

The computation of short range correlation effects is considerably streamlined upon decomposing the nucleon- and 
isobar hole tensors, 



n'^:''^?,") - E ^r\'l,u)L^;iHq,u) + n^^''\q,u)T,Aq,u) 

2 

n( (9, = E ^ {q,u)L^;iHq,u)+ n^^") (g, u) T,, {q, u) 



(8) 



in terms of a complete set of Lorentz structures Llii {q,u) and Tfj_^{q,u). A convenient basis that enjoys projector 
properties was suggested in [l^. We recall the definitions 



(q-u) 



L^^;^\q,u) ^ Ll^^'Hq,u) = q^ 



1 



(q-u) 

5— 91/ - Ui, 

q 

[q-u) 



[q-uf 



■qu -Uy 



(9) 



The presentation of explicit expressions for the longitudinal and transverse nucleon- and isobar-hole loop functions is 
relegated to the Appendix A. The latter follow by simple contraction of the tensors H^y{q^u) with the projectors in 
The results depend on the details of the in-mcdium isobar propagator which will be specified in section 5 and 6. 
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Following [18| we construct the pion self energy in terms of the longitudinal nucleon- and isobar-hole loop functions. 
The self energy can be cast into the form of a sum of 11, 33 and 13, 31 components of an appropriate 4x4 matrix, 



n{q,u) = -47r(l 



) OcffP 

rriN 



11 



13 



33 



31 



where 
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(10) 



(11) 



In (|10p we allow for a background term linear in the nuclear density reflecting a s-wave pion-nucleon interaction. 
Such a term is motivated by the fact that the vertices of ([T]) do not reproduce the empirical s-wave scattering pion- 
nucleon length. At tree-level the vertices H]) lead to a pion-nucleon isospin averaged scattering length of the form 



Anil 



rriN 



-) OttAT 
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rriN 
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1 



2 mA 



(12) 



This leads to o^at ~ —0.09 fm, a significant overestimation of the empirical scattering length of about —0.01 fm 
[23 |. Using the unitarized isobar propagator as implied by the one- loop isobar self energy of section 5 we obtain 
a„jv — +0.00 fm instead, a value significantly reduced and closer to the empirical constraint. In order to correct for 
the remaining slight mismatch we use 6efi — —0.01 fm in pop . 

There are two important technical issues we need to emphasize here. First the application of the longitudinal and 
transverse projectors in ([8|) implies that the loop functions have to satisfy specific constraint conditions. They follow 
from the observation that the polarization tensor Ilf^,y{q,u) is regular, in particular at (7^ = and at (f' — {q ■ u)'^. It 
must hold 



1122(9, ") = nii(g, u) -illi2{q,u) - in2i(g, u) + O [q^) 
n22(9, n) ^ TlT{q, u) + {{q ■ uf - q^) . 



(13) 



The reader may wonder why we discuss this point. After all the integrals ([7]) are finite and the conditions ((T3)) should 
be satisfied automatically. However, we argue in favor of a finite renormalization which is not necessarily compatible 
with ()13|) . A finite renormalization of the isobar-hole loop functions is useful as to suppress the formation of ghosts in 
the pion self energy. The latter may be absorbed into a redefinition of the Migdal's short-range interaction The 
occurrence of ghost causes a severe problem, in particular in a self consistent approach. It implied that the pion self 
energy does not satisfy a Lehman representation anymore. A ghost state is present if the pion self energy has a pole 
for complex energies, i.e. 



L)(w) = det[l-n(^)(w,g)g^] = with 9w7^0. 



(14) 



Note that a function that satisfies a Lehman representation can have poles only on the 2nd or higher Riemann sheets. 
In fact, we observe that such artifacts are avoided typically once a finite renormalization is implemented such that all 
elements Iiij{ijj,q) are bounded for large energies, i.e. 



lim 

; — ^±00 



n,: 



[uj,q)\ < 00 . 



(15) 



As detailed in Appendix A we introduce a finite renormalization for the isobar-hole loop functions by insisting on sub- 
tracted dispersion-integral representations thereof. The construction of the latter was determined by the constraints 
(|13p . We checked that our numerical pion self energies satisfy a once-subtracted dispersion-integral representation to 
reasonable accuracy. In our self consistent simulations we impose such a dispersion-integral representation, where the 
subtraction constant is determined as to find agreement with the direct computation at — = ra,^. 

It is evident that there is a self-consistency issue here. The isobar propagator defining the isobar-hole loop functions 
in ([7]) is a crucial ingredient to evaluate the pion self energy. Since, the isobar self energy depends sensitively on the 
pion propagator itself a self consistent computation is required. The importance of self consistency, as discussed 
above, will be addressed in the second last section when presenting numerical results. 
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IV. ISOBAR PROPAGATOR 

The solution of the Dyson equation ([5]) requires a detailed study of the Lorentz-Dirac structure of the isobar 
propagator. Consider the propagator, Sfj_^{w,u), in the nuclear medium. From covariance we expect a general 
decomposition of the form. 



in terms of invariant functions, S'jj^j (w, u), and a complete set of Dirac-Lorentz tensors u) and Qj^Jj (v, u). For 

latter convenience we introduce 

A suitable basis was constructed in [2^. [25|. enjoying the projector properties 

Qm 9c.p QJJ] = hi Qfy"] , -P^'J = 4i /'[^," . (18) 

This particular basis streamlines the computation of the in-medium part of the isobar self energy significantly. It 
was applied also in In particular the algebra (HH]) illustrates the decoupling of helicity one- half (p-space) and 

three-half modes (q-space). Decomposing the isobar self energy 

S''''(u;,^.) =^ s[^^(t;,«)P[^(z;,«)+^ s[;^^K«)g|:;, («,«), (19) 

into the set of projectors it is straightforward to evaluate the isobar propagator. The Dyson equation (O maps onto 
two simple matrix equations. First, the bare propagator 

needs to be decomposed in terms of the projectors. Second, the six-dimensional matrix E^^"' (w, u) and two-dimensional 
matrix TS'^'^[v,u) have to be evaluated. The final form of the isobar propagator, specified in terms of the invariant 
matrices S''^\v,u) and S"*^'-* (u, w), follows by simple matrix manipulations 

^(P' (w, u) = 5^^' [v, u) \l ~ (w, u) S^"^ (w, m) 1 ~\ 

S^"^ (f , u) = {v, u) \l - {v, u) S^"^ (w, u)j "\ (21) 

The transparent expressions (|2ip rely on the explicit availability of the projector algebra. In order to keep this work 
self contained we review briefiy the set of projectors Pl^^^^{y,u) and Q|^Jj(w,u) introduced in [25j . It is convenient to 
express the latter in terms of appropriate building blocks P±, U±, and L^, i?^ of the form: 



P±{v 




1 ± 


^^2/ 




1 






R^{v 




1 


(f/+(^^, 




u) = 4 


1 

71 


^p(«) ( 



-1 7 • u 



, C/±(«,'u) = P±(«)^==^=_P^(z;), 



(w • u) - Up 



(u • w)^/u^ — 1 



(22) 



For a compilation of useful properties of the building blocks P±, U±, and i?^, we refer to the original work [2f 
The q-space projectors arc 



Q[ii] = ( .9^'' - J ^+ - ^- - P+ P-'^ 
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Q[22] = (.9^" - ^'') - '^^ P+ - ' 

Qfia] = (5''" - ^'^ ^'') ^+ + 5 ^- 

Q[2i] = (.9'"' - ^'^ f/- + i 

+ ^ (l^" P^V + ¥" P+R"") - \L''U-R'' , (23) 



where ^„ = v^/Vv^- It is straightforward to verify (|18p . Using the properties of the building bfocks P±, U±, and 
-^Ai, -R/i [23 reveals that the objects indeed form a projector algebra. 

The p-spacc projectors have similar transparent representations. Following it is convenient to extend the 
p-space algebra including objects with one or no Lorentz index, 



^[11] 
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= c^+, 


^[21] 






^[22] 


= P_ 




^[31] 
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pM 
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= n 




pM 
.[23] 


= P_ 




pM 
M41] 


= y [/_ 


pt^ 


= y^p_ 


pM 

' _[14] 






pM 

.[24] 


= P_ 


V-M , 


pM 
M51] 


= v^' P+ , 


pt^ 
M52] 


= «^ , 


pM 

_[15] 


= n 




pM 

.[25] 


= P_ 
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pM 


= ?)'' C/_ , 


pM 
M62] 


= P_ , 


pM 

M16] 
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pM 

^[26] 


= P_ 




pM 


= P^P+ , 


pt^ 
M72] 


= P^C/+ 


pM 

.[17] 


= n 


P^ , 


pM 
.[27] 


= P_ 


P'^ , 


pA» 
^[81] 


= P^[/_ 


pt^ 

[82] 


= P'' P_ 


pM 

M18] 


= C/4 


P*^ , 


pM 

M28] 


= P_ 


P^ , 



= ^1] - ^] ^2,] ■ (24) 

In the notation of (j24p the indices i, j in ([TBI [T9 l [20 |l run from 3 to 8 in the p-space. The set of identities (fT8|) extends 
naturally 

P[ik] ■ P[ij] = Ski P[ij] , P^.k] Pm = Ski Pi'ij] , -P[^fe] g^v P[ij] (^fe; ^'[y-] , 

Qf,^] 50/3 p£-] = o = P[^] 50/3 g;;]. (25) 

The algebra (j24[) proves convenient in solving various problems. Using the projector formalism we compute the 
in-medium isobar self energy as implied by the interaction vertex ([T]) at the one-loop level in a manifest covariant 
fashion. 



V. ISOBAR SELF ENERGY AND PION-NUCLEON SCATTERING 



It proves convenient to extract the isobar propagator from an appropriately constructed model of the pion-nucleon 
scattering amplitude. Set up in this way all results are induced by expressions already presented in [25| upon the 
application of simple substitution rules. Recall the in-medium Bethe-Salpeter equation, 



T{k, k; w, u) ~ V{k, k; w, u) + 



w ~ l;w,u) = —i S{w — l,u) 



dH 
(2^ 



V(fc, /; w, u) Q{1; w, u) T{1, k; w, u) , 



n{l,u) 



where g, p, p are the initial and final pion and nuclcon 4- momenta and 

w =p + q=p + q, k^^{p-q), k = ^ {p - q) . 



(26) 



(27) 



The two-particle propagator, Q{l;w,u), is specified in terms of the nucleon propagator S{j),u) of ([3]) and the pion 
propagator written in terms of the in-medium self energy n(Z, u) of (|10p . 

In order to generate the isobar self energy S]^^(zi;, u), we introduce the interaction kernel 



P - 

V(fc, k]w,u) = ^ q^ Sl^^iw - u) q^, , 



(28) 
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where we allow for the presence of a vector mean field. The isospin projector is suppressed in p8|) (see e.g. [11]). The 
particular choice (|28p implies a scattering amplitude, which determines the isobar propagator, 5'^i^(w,u), by 



T{k,k;w,u) = -^q^S''''{w,u)q,. (29) 



The system is solved conveniently by decomposing the interaction kernel into a set of projectors, where we apply the 
projectors constructed in terms of the 4-momentum ^ w^ — and rather than and u^: 



^ = E ") 9^ Pm ^) 9. + E ^S] ") 9m 9^ ■ (30) 



For the general case with ^ Yiy the derivation of Vjj^j''''' (w, m) as implied by (^51 ) is somewhat tedious though 
straight forward. The expressions are listed in Appendix B. In the limit — > Sy the expressions simplify with: 

vii) _ T/(p) _ _ /i 1 

^[22] ~ ^[88] 



Mil] 


_ yip) _ I 


/i 1 

ml _ 


y{p) 






3 ml 




^[53] 


- ^[35] ~ ^ 


V3 "4 ™A ' 



"^ir V W + mA 

T^(p) ^ , 2 /a 
[66] ^3 m2 mi ' 

[64] - ^^[46] " V3 m2 mA ' ^ ' 



where only components that are non-zero are specified in (|3ip . A corresponding decomposition is implied for the 
in-medium scattering amplitude 

^ ^Yl P^,-] {v,u)q^ + Y, J^S] (v, u) Qf'^Jj {v, u) q^ , 

T^p) {v,u) = y ) (v, u) \i - j(f ) (w, u) y (w, m) 1 ~\ 

t(«) (w, u) = y (w, u) [i - (w, u) y (u, u)j ~^ (32) 

The scattering amplitude T is determined by the interaction kernel (|3ip and two matrices of loop functions Jj^^^|(w, u) 
and J,^'!i(v,M). Comparing ([32) with (gHl) and (HH) we identify 



J A J A 

> 4 -^Si ■ (33) 

The form of the loop functions can be taken over from [2^, [2^. The evaluation of the real parts of the loop 
functions requires great care. The imaginary parts of the loop functions are unbounded at large energies. Thus power 
divergencies arise if the real parts are evaluated by means of an unsubtracted dispersion-integral ansatz. The task 
is to device a subtraction scheme that avoids kinematical singularities and that eliminates all power divergent terms 
systematically. The latter are unphysical and in a consistent effective field theory approach must be absorbed into 
counter terms. Only the residual strength of the counter terms may be estimated by a naturalness assumption reliably. 
Since we want to neglect such counter terms it is crucial to set up the renormalization in a proper manner. The scheme 
developed in [2^ avoids the occurrence of power-divergent structures and is free of kinematical singularities. 

The loop functions •/[^^j''^ (^^o , ) are expressed in terms of a basis spanned by 13 master loop functions, J„(wo , w ) 
as detailed in [26j . We assume nuclear matter at rest for simplicity. The master loop functions are evaluated by a 
dispersion integral of the form 

Jn{vo,w)^l —1 r J„ (wo,u') , (34) 

J^oc vo - Wo - «e(wo - m) 
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where /x^ 



'AT- 



-kp. We introduce spectral weight functions, AJ„(z;o; vo,w), that depend on 'external' and 'internal' 



energies vq = wq — and Vq. We identify 

AJn{vo;va,w) 



dH 



m% + P 



2(2 7r)3 

K^^K„{l+,vo;vo,w) p^i\v+\,w - I ) 0{+v+) - Q{kF 



(35) 



where the explicit form of the kernels Kn as well as of the counter loops J^{vo,w) are recalled in [2^. The kernels 
are invariant functions of the 4-vectors l^,Vfj_,Vf^ and u^. The spectral density of the pion, pTr{i-u,q), is 



1 



^2 _ ^2 _ „j2 _ ll{uj,q) 
p-^{-^,q) = -^77(^,9) • 



for uj > . 



(36) 



VI. ISOBAR SELF ENERGY IN THE PRESENCE OF VERTEX CORRECTIONS 

The evaluation of the loop functions in the presence of vertex corrections is particularly challenging due to their 
complicated ultraviolet behavior. In Fig. 3 the two types of contributions are depicted graphically in terms of vertex 
functions to be specified below. It is useful to identify a set of master loop functions, in terms of which the full 
loop matrix can be constructed. The latter are rcnormalizcd applying the scheme introduced in the previous section. 
The proper generalization of (|35p is readily worked out. The pion spectral function is distorted by vertex correction 
functions leading to effective spectral densities, which we denote with pah(w, q). For a given spectral distribution we 
introduce 



AJab,n{vo;VO,w) = 

^[k^{1+,vo;vo,w) pab{\v+\,w - I ) 



\v+\ 



a+b 



e(fcf - 1/ 1) 



a+b 



+ 



K„il-,vo;vQ,w) pabi\v^\,w - I ) e(-'y_) | , 



(37) 



where n = 0, ...,12. The kernels Kn{l,va;vo,w) wee identical to those encountered in ((35|) . They are listed in [2^. 
The real part of the loop functions is computed applying the dispersion- integral representation (|34p . A corresponding 
generalization holds for the second term in (|34p . 

We identify the effective spectral distributions, Pab{^, q ) as implied by the diagrams of Fig. 3. The vertex vector 
and tensor may be decomposed into invariants 



T^'{q,u) = q^'T^{q,u) + u^'T^{q,u) , 

r^'^(9, u) = q^ q^'T^^iq, u) + T^2{q, u) + q'' Triiq. u) 

+ r22((?, u) + roo((7, u) , 



(38) 



=*= 



+ 




FIG. 3: Isobar self energy in the presence of short rage correlations. The solid line shows a nucleon propagator in the presence 
of mean fields. The dashed line represents a dressed pion propagator. 



10 



in terms of which wc introduce the spectral distributions 



Lj2 _ g'2 _ ,^j2 _ Il^[uj, q) 



+ ^ab{^^,q) 



(39) 



Applying the techniques introduced in [l^l it is straight forward to express Ti{q,u) and T2{q,u) in terms of the 
longitudinal coupling matrix, .g'^-*, and the loop functions. n(^)(g,w) of (|8I11[) . We obtain: 



Ti{q, u)=\l- f/^'> n(^) {q, u)] .g(-^' H^^) {q, u) 



31 



33 



q ■ u 



^Jq^ - [q- uf 



1 _ g(L) n(i) (g, u)] .g(^) n(-^) (q, u) 



41 



43 



r2(g, u) 



+ ri_gWn(^)(g,^,) 



1 _^g(i)n(^)(g,w) 



41 



43 



(40) 



The matrix Tabiq^ u) probes longitudinal and transverse correlations. As an extension of (jlip we introduce a transverse 
coupling and loop matrix g'-^^ and n'-"^''(g, u). We write 



9ii 9i2 

921 922 



n(T)(^^^)^ /'^^^)(g,^^) 





n(#'^)(g,u) 



(41) 



We derive explicit forms of the tensor vertex 



rii(g,M) = - (x^3^ + , 



(g • uy 



(L) 



q^ — {q ■ u)^ 

ri2(g,w) = r2i(g,u) = - 
r22(g,w) 



X44 



q^ -{q- u) 



v/g2 - (g • u)'- 

(L) (T) 
X44 -X22 



(T) 
2 X22 



q^ — {q ■ uY 

in terms of the longitudinal and transverse correlation functions 

X^^^^\q,u)=\l~g(^-^)n^^^^){q,u) 



. , _ 1 ■ ^ ( (L) _ (T)\ 

^2 _ . „)2 ^X44 X22 ) 

roo(g,w) = X2r ) 



(42) 



(43) 



In the course of deriving the representation ^67\ we made use of the following properties of the correlation functions 

nf{~u:,q) = (-1)'+^- nf{+u:,q) , (44) 
nif (-c.,g') = (-l)^+^n(f (+^,g-), 

r,(-c.,g-) = (-l)''+ir,(+c.,g-), Tab{-io,q) = (-1)^^+^ r„fc(+c., g^) . 

It is left to specify the isobar self energy in terms of the generic loop functions defined by (|37|1 . In a first step a 
matrix of loop functions, J^^'?A, (u, u), is constructed in terms of Jab ni'^^^) detailed in [23|. The latter correspond 



to the projector algebra of section 4. The evaluation of the self energy is analogous to the computation of section 4 
with the slight complication that the effective vertex develops additional structures g^ + gjy, and 5^^. The 
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loops iv,u), which are imphed by the structure 5^ Qi,, contribute hke the previous loops J[^^]'^^(v,m) in das]). The 

implication of the remaining loop functions is readily worked out upon the application of the useful identities 



^ '-MIS] +^26]'' 



V ■ u 



3 V l>2 ^ {^[71] + -^[82] ^ (-^[41] + ^[32])} 

(-^[511 + ^[621 ) ' 



5 ^[11] - + + ^[551 + n'-i 



3 ^[22] - (^[221 + -^^[331 + M661 + MSSl 



'[11] M [55] [77] 
- [22] + M33] 

</ ^12] — V[i2] 3 ^[43] + -^[56] ^ 3 ^[78] 3 >-^[73] ^^[48]'' ' 

9 ^[21] " (<«/[2i] 3 ^[34] + -^[65] + 3 ^[87] ~ "T ^^[Si] + ^[37]> ' 

It is now straight forward to write down the self energies, Yt'^^^^\v , u) . It holds 

2 



ml 

a,b —1 



+ ['^inia] ") + •^2l![ai] 

a=l 
2 



,.6=1 



3 V „2 ^ 
= ^al 5b2 ((5^7 '5j8 - ^i4 (^^3) + (5i5 S^e " ^ (Sir 6j3 + 5i4 Sjsij 

+ 5a2 Sbl (^1 (<5i8 Sj7 - dis dj4) + Ste Sj5 - ^ {Si3 (5j7 + 6is 5j4)^ . (46) 
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VII. NUMERICAL RESULTS AND DISCUSSIONS 



We present and discuss numerical simulations of the pion and isobar spectral distributions at nuclear saturation 
density with the Fermi momentum kp = 270 McV. The results depend on a number of parameters appearing in the 
developed covariant and self consistent approach. These are first of all the scalar and vector mean-field shifts of the 
delta, and S^, as well as the Migdal parameters g'n^ g'12 and (722- One should also consider medium induced 
changes in the coupling constants /a and /^r, although it is usually assumed that for nuclear densities they do not 
depart significantly from their vacuum values. The nucleon mean-field parameters 'Sfj and Sy, which model nuclear 
saturation and binding effects, are al so p rone to variations in different models. We use the values = 350 MeV and 
Y.^ = 290 MeV also assumed in Hi, [23|. 

As a guide we consider the values of the above parameters used in earlier computations, but having in mind their 
scheme dependence. We focus on variations around a parameter set which has been shown to lead to a reproduction 
of the nuclear photoabsorption cross-section in the delta excitation region [23] . The latter study builds on the 
self-consistent approach developed in this work. It is the first work that considers photabsorption in the presence 
of short-range correlation effects in the jinr, j N A, -firNA, ttTVA and tt N N vertices. Electromagnetic gauge 
invariance is kept as a consequence of a series of Ward identities obeyed in the computation. In particular the 
interference of the in-medium s-channel isobar exchange and the t-channel in-medium pion exchange is considered. 
We refer to the details of that work which provides the following parameter set: 

S]t = -0.25GeV, S^^-O.llGeV, g'n^l.O, g'12 = 0.4, g^2-0.4, (47) 

together with an in-medium reduction of the /a coupling by 15% but an unchanged value for f^. According to 
the in- medium reduction of the /n coupling is quite small (less than 6% at saturation density). This is in line with 
our finding that the photoabsorption data does not require any in-medium change of /at. The values of the Migdal 
parameters in (|47p is within range of the various sets used in the literature. Though in the recent work by Hees and 
Rapp [2^ large values of g'12 and g'22 are excluded in their non-relativistic scheme, this is not the case in our more 
microscopic and relativistic approach. Values for 17^2 and g'22 as large as in ()47|) generate a width for the isobar in f23l | 
that would be incompatible with the photoabsorption data as computed in [23| . However, it is reasonable to expect 
that such a condition is altered by a possible in- medium reduction of /a- 

We recall from 0,123] that the actual position of the photoabsorption peak is a subtle effect of short-range correlation 
effects and the in-medium isobar properties. The peak of the isobar spectral distribution does not translate directly 
into the maximum of the photoabsorption cross section. The pion and isobar properties as implied by (|47l) are shown 
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FIG. 5: Contour 
plots of the pion 
spectral function us- 
ing the parameters 
(gTJ. The upper 
left figure gives our 
full result, the upper 
right figure follows 
if the vertex correc- 
tion in the delta self- 
energy are neglected. 
The lower left figure 
follows if the free- 
space isobar propa- 
g gator is used in pup . 
The lower right figure 
shows the impact of 
reducing the nucleon 
mean fields. 
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in Fig. [4] for nuclear saturation densities by solid lines: at zero momentum the isobar receives an attractive mass 
shift of about 50 MeV. A value amazinglyclose to the range obtained in 01 but in stark contrast to the small and 
repulsive mass shift obtained recently in [23|. For the isobar we restrict the discussion to the two main components 
because they dominate the resonance region. Please note however, that the proper inclusion of all other components is 
essential to ensure the cancelation of kinematical singularities on the light-cone. We observe a significant splitting of 
the p- and q-spacc modes at nonzero momentum. The medium effects are stronger for the q-space (hclicity 3/2) than 
they are for the p-space (helicity 1/2), where we obtain a less pronounced broadening and smaller shift in the position 
of the peak at larger momentum. Note that the nuclear photoabsorption data probe dominantly the helicity 3/2 
mode. These finding are in qualitative agreement with the results of [7| that were based on a perturbative and non- 
relativistic many-body approach. It should be pointed out, however, that the pion spectral function corresponding to 
the approach of Q differs decisively from the one predicted by our approach. Though a direct comparison is difficult, 
since Oset et al did not provide figures for the pion spectral function, an indirect comparison may be possible. We take 
the more recent work of Ramos and Oset [l^ , which provides explicit results for the pion spectral distribution. The 
strength in the soft pion modes as shown in Fig. [3] is much suppressed as compared to an in-medium pion considered 
realistic in [3^ . Also a comparison of our pion spectral function in Fig. [3] with other recent results [l^, [13, HI] show 
significant and systematic differences at small and intermediate momenta. 

Before we discuss a variation of parameters around the central values (j47[) wc examine the effect of various ap- 
proximations all based on the parameter se t (1471). First we consider the effect of neglecting short-range correlation 
effects in the isobar self energy as done in |12|. Il9l. [20| . In Figs. S] and [5] the quality of such an approximation is 
scrutinized. Though the Migdal parameters enter in a decisive manner in the computation of the pion-self energy via 
(1 1011, th e isobar properties are a functional of the pion self energy only. As studied in great detail in the previous works 
[l2l. fiol. [20|. the self consistent treatment of the pion and isobar properties is an important and significant effect even 
in the absence of vertex corrections for the isobar. The upper left panel of Fig. [5] shows the contour lines of the pion 
spectral function as obtained in the fully self consistent computation. If one neglects the ttNA vertex correction in 
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the isobar self-energy as discussed above, the contour hues in the upper right-hand panel arise. A more quantitative 
illustration is offered by a comparison of the solid and dashed lines in Fig. [H The Figs. 2] and O document the 
importance of the vertex corrections in the isobar self energy. Most significant are the effects on the isobar spectral 
distribution as shown by the solid and dashed lines of Fig. ID The consistent consideration of short range correlation 
effects leads to a significant attractive mass shift and a reduction of the width for the isobar. It is interesting to 
observe that it appears well justified to treat the vertex contributions in the isobar self energy in perturbation theory. 
We find that the evaluation of the vertex bubbles of Fig. [3] with a free-space isobar propagator leads to results that 
can barely be discriminated from our full results. Recall, however, that a corresponding attempt for the short-range 
bubbles in the pion self energy would fail miserably. This is illustrated by the lower left-hand panel of Fig. O where 
the pion spectral function is shown as it is implied by the free-space isobar together with the Migdal parameters 
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of (j47p and our in- medium value for /a- In particular the width of the low- momentum main pion mode would be 
underestimated. 

We now turn to a variation of the parameter set. In the lower right-hand panel Fig. [5] the effect of using smaller 
scalar and vector mean fields for the nucleons is illustrated. The contour lines were obtained with T,g = 175 MeV and 
Sy = 115 MeV. In Fig. [S]the variation of the value chosen for the ttNA coupling constant /a is investigated. The 
reason for considering the departure from the vacuum value is that a detailed study j27j of nuclear photoabsorption 
strongly favors such a change, more precisely a reduction of the /a coupling by about (10 — 15)% at nucleon densities 
close to saturation. As expected a reduction of the coupling by leads to a reduction of the isobar width. In the pion 
spectral function we can, at least for intermediate momenta, distinguish three branches. These are the main pion 
mode as well as the particle-hole and A-hole excitation. At about 0.3 GeV momentum we observe the level crossing 
between the main pion mode and the isobar-hole excitation. Decreasing /a reduces the strength of the isobar-hole 
branch and in addition due to the narrower isobar that mode becomes better visible. 

Next we study the influence of Migdal's g'22 parameter. Varying its value from 0.2 to 0.5 we arrive at the results 
shown in Fig. I VIII The effect of changing g'22 is subtle since it influences the dressing of the isobar through the ttNA 
vertex correction and also the pion self energy by affecting the isobar-hole loop contribution. Increasing the value of g'22 
softens the isobar and decreases its width, which compensates in part the reduction of the isobar-hole-loop contribution 
to the pion self energy. All together the resulting change in the pion spectral function is modest. We note that a 
variation of g'12 is quite similar to that of 322- ^ variation of g'n just affects the nucleon-hole contribution. Lowering 
g'li makes the nucleon-hole branch of the pion larger, which in turn somewhat increases the isobar broadening. 

We conclude with a discussion of the influence of the isobar mean field parameters. Results are shown for two 
parameter sets, which induce the same energy shift at zero momentum. Next to our standard set we use a set whose 
scalar mean field is put to zero and the vector part provides the net repulsion of 0.14 GeV at zero momentum as is 
implied also by (|T7)) . The effects can be found in Fig. 8. Without the scalar mean field we obtain less attraction at 
nonzero momentum and in addition the width of the isobar is significantly increased at larger momenta. This implies 
a smaller contribution of the isobar-hole state to the self energy of the pion as shown in the lower right-hand panel of 
Fig. 8. 



VIII. SUMMARY 



A detailed study of pion and isobar properties in cold nuclear was presented. A fully relativistic and self-consistent 
many-body approach was developed that is applicable in the presence of Migdal's short range correlations effects. 
Nuclear saturation and binding effects were modeled by scalar and vector mean fields for the nucleon. The novel 
subtraction scheme, that was constructed recently by two of the authors and that avoids the occurrence of kinematical 
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singularities, was used. Unlike in previous studies no soft form factors for the ttTVA vertex were needed. For the first 
time the ttNA vertex corrections as dictated by Migdal's short-range interactions were considered in a relativistic 
and self consistent many-body approach. The latter were found to affect the isobar and pion properties dramatically. 
Using realistic parameters sets we predict a downward shift of about 50 MeV for the A resonance at nuclear saturation 
density. The pionic soft modes are much less pronounced than in previous studies. 

Further studies are needed to consolidate our results. In particular an application to the pion- nucleus problem and 
the pionic atom data set would be useful to further constrain the parameter set. Our computation may be generalized 
to study effects of finite temperature. 
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Appendix A 

We present explicit representations for the nucleon- and isobar-hole loops introduced in ([7]) for the case of nuclear 
matter at rest with li^ — (1,0). The longitudinal and transverse nucleon-hole loop functions are: 



Yr(Nh) , _ In_T> f 



+ (48) 



where = (t^,g), Po = \J mj^ + p^ + Sy and 



fC^Nh) _ 2 r^{Nh) _ (Nh) _ 

^11 — ^"ZjY, J\i2 —^21 



^.(^/.) ^ |^2p-2 + (po - Sy) + p • gj + 2 mj 

= 3 < + c (po - Sy) - P- q- \ (kT' + kT^) ■ (49) 

For a bare isobar propagator, Sq'^{w) as given in ([5|), the longitudinal isobar-hole loop functions were computed 
already in [l^. We present here longitudinal as well as the transverse loop functions: 



.(A.). ^ ifl f"^ d^p Shj^'^ {m^m^ + ml + ip-q)) 

^'^''^^ QmlJo 2 (po - Sy) (2^)3 2p ■ q + q^ - + m], + i e 



+ {-iy+^{q,^-q,), 
r^iAh) _ . _ {q^ +P-qf j^(Ah) _ , juJ \p\ cos{q,p) - \q\poy 

-f'^12 --^21 -'>'— — —[\<l\Po~'^\p\cos[q,p)) 
q 

i^r = 2-^4(Am^^f'), (50) 

where = {i-u,q), Po = \/ mj^ + p^ + Sy. Both representations (|48[ [50)) are compatible with (fT5|) . On the other 
hand, only ([48]) is consistent with p^ . The asymptotic behavior of the isobar- hole loop as given in (|50p is at odds 
with the condition 

To derive the general results for the isobar-hole loop functions it is advantageous to choose a representation slightly 
different to (|5D|) . We write 

nil {^,1) = — ni (w,(7) , 



q- 

('^:5') = — F====== — n 'w,g-n^ '{uj,q) 

^q^ - [q ■ u)2 V 9^ 



^22 n^.g) = g2„ (g.^)2 I (w,g)-2n^ ^(w,g) 



q-u / 



n<^^''^(-,gn = i (nf'^ (-,9--)- nir (-,<?) -nr '(-,?)) • (5i) 



The merit of the representation ([?T|) lies in its simple realization of the constraint equations (fT51) . The first condition 
is satisfied for any functions Tli{LO,q) that are regular at q^ = 0. The second equation in (|13|) implies the following 
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constraint, 



1 



1 



U3{u;,0) = — ni(a;,0), n2(w,0) ^-Ui{uj,0), 

Ui{u;,0) = 3U22ii^,0) + Uii{u;,0), (52) 

where we boosted into the rest frame of nuclear matter for convenience. Based on the representation p6p we define 



n 



{Ah] 



fi 



3 ml Jo 2(po-Sy)(27r)3 7_ 



+ 00 



duj (U!\^i sign (w) 35*1'^''^ (w, 



TT V(D/ 



tjj — tjj — iLO e 



(53) 



where po = +?'^ + a-nd £1,3,4 = and (.2 = 1- Furthermore ni,4 = 2 but 712 = 1 and 713 = 0. We assure 

that the definition ([53|l leads to a polarization tensor compatible with all constraints [T5|l . This is a consequence 
of specific identities the integral kernels enjoy (see [56]). 

The integral kernels, S^^^'^\q,p,u), required in (|53p are covariant functions of the 4-momenta qf_t,Pf_t and m^. Their 
evaluation requires the contraction of the isobar propagator, S^y{p-\-q, u), with the q^ and (see (I3[9])). We express 
the 4- vector u^, in terms of and Xp(w,u), 



= -\J {v ■ uY - lX^{v,u) + {v ■ u)Vf^, 



(54) 



since the contraction of the isobar propagator with and Xfj^{v,u) leads to more transparent expressions. In 
particular we can take over the results from (26j , where contractions of the isobar propagator with the latter 4- vectors 
were computed already. The results were decomposed into the extended algebra of projectors p3l [24]) introducing 
the invariant expansion coefficients S^°'j-^{v,u) and S^^^'^\v,u) with a,b = v,x. 



We present the integral kernels of (|53p. which have transparent representations in terms of the invariant functions 

M 



introduced in ^W\i and ci^'?\q;v,u) of (26|. We establish 



c(A/i) _ M qip) , Jq) oil) 

~ "^fo] ^ "^fo] ' 



(A/i) _ (p) 
2^ 2^ "^[y] 

i=l j=3 
2 



-t(a;a;) 



-(^.u)V(«-)^-i(^g]^+^f;;f 
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(P) c(9) 
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(55) 



A straight forward computation reveals that the kernels S^'^^^ are correlated at vanishing 3-momcntum q =0. In 
this case it holds 



S. 



(Ah) 
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(Ah) 



s. 
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MAh) 



(56) 



(Ah) 



^ dq^ 



g=0 



where we assumed an angle average, i.e. the presence of dflff. 
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